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ABSTRACT 

High-frequency quasi-periodic oscillations (QPOs) appear in general-relativistic magnetohydrodynamic sim- 
ulations of magnetically choked accretion flows around rapidly rotating black holes (BHs). We perform polar- 
ized radiative transfer calculations with our ASTRORAY code in order to explore the manifestations of these 
QPOs for Sgr A*. We construct a simulation-based model of a radiatively inefficient accretion flow and find 
model parameters by fitting the mean observed polarized source spectrum. The simulated QPOs have a total 
sub-mm flux amplitude of under 5% and a linearly polarized flux amplitude of up to 2%. The oscillation period 
T w lOOM w 35 min corresponds to the rotation period of the BH magnetosphere that produces a trailing spiral 
in resolved disk images. The total flux signal is statistically significant over noise for all tested frequencies 
87 GHz, 230 GHz, and 857 GHz and incHnation angles 10°, 37°, and 80°. The non-detection in the 230 GHz 
Sub-Millimeter Array light curve is consistent with a low signal level and low sampling rate. The possible 
presence of such magnetospheric QPOs in Sgr A* will be better tested with Atacama Large Millimeter Array. 
Subject headings: accretion, accretion disks — black hole physics — Galaxy: center — instabilities — mag- 
netohydrodynamic s (MHD) — radiative transfer 



1. INTRODUCTION 

Quasi-periodic oscillations (QPOs) in the emission from 
black hole (BH) accretion disks and jets are found in sys- 
tems with both stellar mass BHs ( Remillard & McClintocH 
2006") and supermassive BHs (SMBHs) ( Gierlihski et aT| 
2008; Reis et al. 2012). High-frequency QPOs (HFQPOs) 
with period T that is similar to the orbital period at the in- 
nermost stable circular orbit (ISCO) potentially probe the re- 
gion close to the BH event horizon, offering a chance to test 
accretion and jet theories in the strong gravity regime. There 
have been multiple claims of HFQPOs from the SMBH Sgr 
A* in the Milky Way center, which then provides a unique 
opportunity to study HFQPOs up-close. Henceforth, we set 
the speed of light and gravitational constant to unity (c = 1 
and G = 1), such tha t IM = 21 s for Sgr A* with BH mass 
Mbh = 4.3 X lO^Mo (iGhez et al.ll200l iGillessene t al. 2009). 

A HFQPO period commonly claimed to be seen in Sgr 
A* is r = 17min = 4 8M. This was id entified in the K 
ba nd during a flare (iGenzel et al.l 120031) . in 7 mm data 
by lYusef-Zadeh et al.] (120111) . and in images obtained with 
very l ong baseline interferometry (VLBI) by Mivoshi et id] 
(l201lh . Periods of T = 28min = 79M (iGenzele t al. 2003^ 
T = 23min = 65M, and T = 45 min = 127M iTrippe et al. 
120071; iHamaus et al. 120091) were rep orted in IR observations 
of flares. However. IDo et al.l (l2009f) analyzed a long K band 
light curve and did not find power spectrum density excess to 
be statistically significant at any period. A statistically signif- 
icant longer period o f T = 2.5- 3hrs « 470M was repo r ted by 
iMauerhan et al.1 (2005) in 3 mm data. Mivoshi etal] ( 120111) 
found a range of periods from 17 min to 56 min claiming to 
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see the spirals with multiple arms, but did not perform signif- 
icance analysis. 

The physical origins of HFQPOs are highly debated. 
Models are often based on the ISCO orbital frequency, 
the epicyclic frequency, and the fr equencies o f vari- 
ous pressure and gravit y mod es (see iKatol 1200 IL l2004t 
iRemillard & McChnt ock '2006'; 'Wagoneij |2"008rfor an 
overview). Some underlying physical origins may be 
beat frequencies (fvan der Klis 2000), resonances between 
flow modes and normal frequenci es in general relativ- 
ity (GR) (Abramowicz & Kluzniala 1200 Ih . trapped oscil- 
lations (Nowak & Wagoneri il991i) . parametric resonance 
(Abramowicz et al. 2003), and disk magnetospheric oscilla- 
tions (Li & Narayan 2004). 

Such analytic methods have been followed by searches 
for QPOs in magnetohydrodyna mic (MHD) simu lations. 
Non-GR 2D MHD simulations by iTagger & Melial (|2006); 
iFalanga et al.l (120071) exhibited spiral patterns of Rossby 
waves detectable in simulated light curves. Non- GR 3D 
MHD simulations of thick accretion disks by Cha n et al] 
(I2009h developed QPOs with a period T = 39M in sim- 
ulated X-ray light curve. A s imilar search in thin disk 
simulations (Reynolds & Mille] l2009h yielded no QPOs. 
The 3D GRMHD s imulations and radiative transfer by 
ISchnittman et al.l (12006 ) revealed weak transient QPOs. Sim- 
ilar simulations by iDolence et aT (2012) post-processed with 
radiative transfer code grmonty (Dolence et al. 2009) showed 
a spiral structure producing oscillations with periods T = 
6-9min= 17-25Min simulated NIR andX-ray light curves. 
GRMHD simulations of tilted disks produce low quality os- 
cillations in dynamical quantities with T k, \1QM = 1 hr 
dHenisev et al. 2012). 

Such MHD simulations start with a weak magnetic field, 
which is amplified by the magneto-rotational instability 
(MRI) that generates incoherent turbulence. However, when 
magnetized gas falls onto a BH, the disk can become sat- 
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urat ed with more magnet ic flux than the MRI can gener- 
ate jMcKinnev et alj|2012l) . 3D GRMHD simulations of ra- 
diatively inefficient accretion flows (RIAFs, as appUcable 
to SgrA*) with ordere d magnetic flux were performed by 
iMcKinney et alJ ( l2012h . The resultant magnetically chocked 
accretion flow (MCAF) has a BH magnetosphere that signif- 
icantly affects the sub-Keplerian equatorial inflow. The sim- 
ulations showed high-quality disk-BH magnetospheric QPOs 
in dynamical quantities with an m = 1 (one-arm) toroidal mode 
and a rotating inflow pattern in the equatorial plane. Such 
magnetospheric QPO s are robust to signific ant tilt between 
the disk and BH spin dMcKinnev et alj|2013h . 

In present paper, we quantify the QPO signal and its sta- 
tistical significance in simulated Sgr A* light curves for flow 
models based on GRMHD simulations of MCAFs. In Sec- 
tion|2]we describe the 3D GRMHD simulations and the appli- 
cation to Sgr A*. We perform GR polarized radiative transfer 
calculations with our ASTRORAY code and fit Sgr A* mean 
polarized spectrum. We find the best-fitting accretion rate M, 
electron temperature Te, and the BH spin inclination angle 0, 
which is coincident with the viewing angle. In Section [3] we 
describe timing analysis. We study the light curves of the 
best-fitting model viewed at different 9. We find QPOs to 
be statistically significant in total and linearly polarized (LP) 
fluxes. While the largest absolute flux variations are achieved 
in total flux, the largest relative variations are exhibited in 
LP flux. We image a correspondent equatorial plane spiral 
wave. In Section |4] we compare the simulated QPO ampli- 
tudes and periods with those of observed Sgr A* oscillations. 
We argue that the non-detection in 230 GHz Submillimeter 
Array (SMA) data is consistent with the low level of simulated 
QPOs, but future Atacama Large Millimeter Array (ALMA) 
observations may detect sub-mm Sgr A* oscillations. 

2. SGR A* MODEL BASED ON GRMHD SIMULATIONS 

2.L GRMHD Simulations 

The initial gas reservoir is a hydrostatic torus 
dGammie et al.l l2003h within which magnetic field loops 
are inserted. The MRI action on the initial field leads 
to MHD turbulent accretion that eventual ly causes mag- 
netic flux to saturate near the BH jMcKin nev et al. 2012). 
We focus on a simulation with a dimensionless spin 
a, = 0.9375, which is similar to a* k, 0.9 favored in 
simulation-based m odeling of Sgr A* spect r um and the 
emitting region size ("Moscibr odzka et al]|2009l: iDexter et al.l 
120101: Shcherbakov et al. 2013)^ The simulation is per- 
formed in spherical coordinates (r,9,(j)) with resolution 
Nr X Ng X Ncj, = 272 X 128 x 256. It reached a quasi-steady 
state by time t = 8,000M and ran tiU f = 28,000M. In steady 
state near the BH event horizon, the sub-Keplerian inflow is 
balanced against the BH magnetosphere resulting in vertical 
compression of the disk. 

The BH magnetosphere and disk exhibit QPOs in dynam- 
ical quantities such as the magnetic field energy density. A 
toroidal wobbling mode with m = 1 is eminent in the jet po- 
lar region and disk plane. They were identified with pattern 
rotation of the BH magnetospheric region pierced by the in- 
falling matter streams. The strea ms form due to magneti c 
Rayleigh-Taylor instabilities (e.g. IStone & Gardiner! 120071) . 
The pattern rotates with an angular frequency fif « 0.2J7//, 
where Hh = fl*/(2r/f) is the BH angular frequency and rn = 
(I + ^yl-al)M is the horizon radius. The angular frequency 
flp is similar to the rotation frequency ~ 0.270// of the field 



lines attached to the BH a t the equatorial plane in par aboloidal 
magnetospheric solution dBlandford & Znaiefll 19771) . 

2.2. Sgr A * Accretion Flow Model 

We use this 3D GRMHD simulation of an MCAF to model 
the Sgr A* accretion flow. We follow [Shcherbakov et all 
(2012) to define the electron temperature and continuation 
of quantities to radii r outside of 50M. A power-law ex- 
tension of density to r > 50M is n c>c r~^, while the proton 
temperature is continued as Tp oc r"' . The magnetic field is 
extended as fe ex JnTp oc r'"'"*^'/^ to preserve a constant ra- 
tio of the magnetic field energy density to the thermal en- 
ergy density. The slope /3 is found by connecting the known 
density at r = 3 x lO^M to the density in the inner region 
(Shcherbakov & Baganoff 2010). Correct simultaneous evo- 
lution of the simulations and the radiation field is considered, 
despite radiative transfer is conducted in post-processing. 

We focus on the accretion disk as the source of emission 
from Sgr A* and will consider jet emission (e.g.. lFalcke et alj 
120041) in future studies. The simulation's matter densities are 
artificial near the polar axis because matter is injected in or- 
der to avoid an exceedingly high local ratio of magnetic en- 
ergy to rest-mass energy that is difficult for GRMHD codes to 
evolve. The injected material does not change flow dynamics 
because it's locally energetically negligible. However, a small 
amount of hot matter in t he polar r egion can shine brightly 
as revealed by Moscibrod zka et all (l2009h . The matter den- 
sities are zeroed out in a bipolar cone with an opening angle 
B = 26°. If that artificial matter was not removed, then none 
of our models w ould be consistent wit h the observed image 
size at 230 GHz dDoeleman et alj|2008l) and the observed po- 
larized Sgr A* spectrum. 

Radiative transfer is performed with our ASTRORAY code 
dShcherbakov & Huang 2011; Shcherbakov et al. 2012). We 
compute radiation over a quasi-steady simulation period be- 
tween t = 8,000M and t = 28,000M. Following our previous 
work, we fit the total flux of Sgr A* at 87 - 857 GHz, the 
LP fraction at 87 GHz, 230 GHz, and 345 GHz, and the cir- 
cular polarization (CP) fraction at 230 GHz and 345 GHz. 
We vary the heating constant C, which determines the elec- 
tron temperature Tg close to the BH, the accretion rate M, 
and the inclination angle 9. Comparing the mean Sgr A* 
spectrum with the mean simulated spectra we reach x^/dof = 
1 .55 for dof = 9, which is a better fit th an our prior work 
based on weakly magnet ized simulations dPenna et al.ll2010l : 
[Shcherbakov et al] 120121) . The correspondent values of pa- 
rameters are J^ = 3.2 x 10^" K at 6M distance from the center, 
M = 1 .0 X 10"^Moyr"\ and 9 = 31°. We then perform timing 
analysis of the model, which reproduces Sgr A* spectrum. 

3. TIMING ANALYSIS 

3.L Oscillations in Light Curves and Images 

Let us first demonstrate oscillations in the light curves. In 
Figure[T]we show the light curves between times t = 25,500M 
and f = 26, lOOM for the best-fitting model with the inclination 
angle 9 = 31°. The light curves are shown at three frequencies 
with different optical depth: radiation at 87 GHz is optically 
thick, the optical depth at 230 GHz is about t ^ I, while ra- 
diation is optically thin at 857 GHz. The total flux (top panel) 
shows regular oscillations with the amplitude AF « 0.05 Jy 
at 87 GHz and AF w 0.15 Jy at 230 GHz. Fluctuations at 
857 GHz with the amplitude AF « 0.2 Jy are less regular. 
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Fig. 1. — Fragments of the light curves for the best-fitting model with 
the inclination angle d — 'il° at optically thick 87 GHz (blue solid line), at 
230 GHz (red dashed line) with the optical depth about unity, and at optically 
thin 857 GHz (green dotted line). 



The LP fraction fluctuates at 2% level at all three frequencies, 
which translates into relative variations of up to 50% and ab- 
solute LP flux variations of AF w 0.06 Jy. The LP and CP 
fractions and the electric vector position angle (EVPA) ex- 
hibit substantial variations over long timescales at 87 GHz. 
Variations of the EVPA at 230 GHz and 857 GHz are about 
5° - 10°. The CP fraction osciUates by 0.3% at 87 GHz and 
by 0.15% at higher frequencies. Absolute CP flux variations 
are AF w 0.005 Jy at 230 GHz. 

In Figure |2] we show how the amplitude of oscillations de- 
pends on the inclination angle 6 at 230 GHz. Shown are the 
light curves for the best-fitting inclination angle 9 = 37° (solid 
line), for almost face-on 6* = 10° (dashed line), and for almost 
edge-on 9 = 80° (dotted line). The total flux shows the same 
amplitude AF ?» 0.15 Jy regardless of 9. The edge-on and 
the best-fitting cases produce comparable variations of the LP 
fraction, while cancelations of the polarized fluxes emitted in 
different regions lower both the mean and the fluctuation am- 
plitude of the LP fraction in the face-on case. Correspond- 
ingly, the EVPA fluctuates dramatically in the face-on case. 
The CP fraction oscillates at 0.5% level in the edge-on case, 
while the other cases exhibit 0.15% amplitude of CP oscilla- 
tions. 

The face-on accretion flow images with = 10° are shown 



Fig. 2. — Fragments of the light curves at 230 GHz for the best-fitting Te 
and M and several inclination angles 6: best-fitting 9 = 37° (solid), face-on 
e = 10° (dashed), and edge-on 6» = 80° (dotted). 

in Figure [3] We put the time series of the total intensity 
images in the top row and the time series of the LP inten- 
sity images in the bottom row. The total intensity images 
show a clear one-arm spiral, which rotates with a period of 
T « lOOM w 35 min. A spiral of LP intensity is offset from 
the spiral of total intensity, as the region of the brightest to- 
tal intensity exhibits the strongest LP cancelations. Our total 
intensity spiral looks similar to that in Dolence et al.l (1201 2|) . 
despite different angul ar velocities. Note th at we report the 
intensity images, while iDolence et aTl (120121) showed the im- 
ages of the dynamical quantities. 

3.2. Statistical Analysis 
Let us quantify the signifi cance of QPOs. Following 



iPapadakis & Lawrence! (Il993h we start with an autocovari- 
ance 
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Fig. 3. — Images of a face-on disk at 230 GHz: total intensity images (top row) and LP intensity images (bottom row). Strokes indicate the EVPA direction. 
The rotating spiral pattern is clearly visible. 



In our analysis At = AM, which appears large enough to 
avoid aliasing at periods T > 50M. All periodograms are 
log-smoothed to O.OSdex following iPapadakis & Lawrence! 
(fT993i) . 

Determination of the statistical significance of QPOs in- 
volves comparison of the simulated periodogram with ran- 
dom noise periodogram s. We follow the procedure in 
iTimmer & Koeni j (Il995h to generate random noise. We em- 
ploy a log-smoothed to 3.0dex periodogram of the simulated 
light curve to normalize the noise power, which eliminates 
power spectrum peaks and introduces less bias compared to 
fitting the simulated periodogram with a powerlaw. We draw 
the Fourier transform of random noise from a normalized 
Gaussian distribution and perform the inverse Fourier trans- 
form to generate the light curve. We then compute random 
noise autoco variance and periodogram. We find a 3cr signifi- 
cance curve based on 2,592 random noise samples. 

Periodograms with correspondent 3 a significance curves 
are depicted in Figure |4] for a number of frequencies and in- 
clination angles. The top curves in each panel are for al- 
most face-on inclination 6 = 10°, the middle curves are for 
9 = 31°, and the bottom curves are for almost edge-on 9 = 80°. 
QPOs with a given period are statistically significant, if a peri- 
odogram rises above the 3fT curve at that period. The total flux 
and the LP fraction periodograms exhibit statistically signifi- 
cant peaks with the period T « lOOM at most studied frequen- 
cies and inclination angles. However, a face-on disk generally 
shows weaker QPOs, which are not significant at 857 GHz. 
LP fraction oscillations are weak for a face-on disk at all fre- 
quencies due to random cancelations of the LP. The total flux 
and the LP fraction at 230 GHz show the strongest signal in- 
dicating that current active efforts to observe at 1 .3 mm wave- 
length are the most promising for QPOs. Oscillations with a 
period T = lOOOM « 4 hr are detected with 3cr significance 
in LP fraction. We characterize presence of oscillations and 
stability of an oscillation period by a spectrogram in Figure|5] 
The spectrogram indicates that most of the time oscillations 
with a period T = 90- IQQM are present. However, no os- 
cillations occur around time t = 22,000M, when a burst of 
accretion is caused by weaker magnetic field. 



4. DISCUSSION AND CONCLUSIONS 
4.1. Summary and Comparison to Previous Work 

In this paper we report QPOs in simulated Sgr A* light 
curves for a model based on the state-of-the-art 3D GRMHD 
simulations of magnetically chocked RIAFs. By adjusting the 
parameters we find a fit with x^/dof = L55 for dof = 9 to 
the mean sub-mm source spectrum. Correspondent simulated 
total flux light curve shows regular oscillations with the pe- 
riod T w lOOM « 35 min and the ampHtude AF w 0.15 Jy at 
230 GHz. Less regular fluctuations with AF « 0.2 Jy are seen 
at 857 GHz. Weaker oscillations with AF « 0.05 Jy are seen 
at 87 GHz, which probes the optically thick emission from 
^ lOM radius. The LP fraction exhibits periodic modulations 
at 50% relative level, but the absolute LP flux amplitude is 
only about AF « 0.06 Jy. QPOs are significant at a 3(7 level 
in the total flux light curves for all tested inclination angles 
10°, 37°, and 80° and frequencies 87 GHz, 230 GHz, and 
857 GHz, while the LP fraction shows less prominent QPOs 
at 87 GHz and in a face-on case. 

Our main T « 35 min period is longer than the claimed ob- 
served 17-20 min pe riod of Sgr A* , while the simulated pe- 
riods T = 6-9 min in Dolence et alj (120 12h are shorter Their 
simulation has the same BH spin a* = 0.9375 as in our simu- 
lation, but Dole nce et alj (|2012) simulations evolve to below 
saturation levels for the magnetic flux near the BH and to a 
somewhat thin disk with height-to-radius ratio of H/R ^ 0.2. 
The resultant MRI-dominated accretion flow has Keplerian 
rotation and shorter QPO periods. They did not identify their 
QPO mechanism, although they noted their turbulence is un- 
resolved (iShiokawa et al. 2012) a nd this might lead to arti- 
ficial QPOs dHenisev et alj 120091) . Our MCAF model has 
sub-Keplerian rotation and QPOs drive n by the rotating BH 
magnetosphere's actions on the disk dLi & NaravanI l2004b 
that leads to longer periods. Our relatively thick disk with 
H/R ^ 0.6 is similar to expected for a RIAF, and ou r sim- 
ulations well-resolve disk turbulence ( McKinnev et al.ll2012h 
suggesting the QPOs are robust. Based upon these works, Sgr 
A* QPOs might be explained by a* = 0.9375 with an interme- 
diate gas rotation rate, magnetization, or H/R. 

Flow cooling, whose marginal importance for Sgr A* was 
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Fig. 4. — Periodograms of the total flux and the LP fraction Hght curves for the frequencies 87 GHz, 230 GHz, and 857 GHz and the inclination angles 9 = 10° 
(top curves in each panel), 9 = 37° (middles curves), and 6 = 80° (bottom curves), while Te and M are fixed at their best-fitting values: simulated periodograms 
(dark black solid lines), 3(t significance curves (Hght red solid lines), log-smoothed to 3.0dex simulated periodograms (dark black dashed lines), and geometric 
means of random noise periodograms (light green dashed lines). 



suggested by iDrappeau et al.l (120 12h . can self-consistently 
choose the disk thickness H/R in simulations. In MCAFs, the 
steady-state BH horizon magnetic flux has a positive corre- 
lation with H/R (iMcKinney et alJ[2012h . so cooling can lead 
to less sub-Keplerian rotation rates and a weaker magneto- 
sphere, and then the the QPO period from MCAFs could be 
comparable to claimed for Sgr A*. 

4.2. Observing QPOs in Sgr A* 

We showed that QPOs, though significant, have maximum 
amplitude of 5% or AF ~ 0.15 Jy in the sub-mm. Low 
sensitivity and low sampling rate of current sub-mm instru- 
ments might prohibit o bservational detection of such oscil- 
lations dMarroneJ 120061) . SMA achieves an accuracy of 5% 
and sampling every 10 min at 1.3 mm with a correspon- 
dent 20 min Nyquist period (iMarrone et al.ll2008l) . Thus, a 
weak signal with T ^ 30 min period can be readily masked 
by ahasing and noise in SMA data. Caltech Submillimeter 



Observatory (CSO) is comparable with SMA in sensitivity 
and sampling rate (Marr one et al.ll2008l) . Very Large Array 
(VLA) provides exceptional signal-to-noise ratio and sam- 
phng rate to detec t QPOs at higher wavelengths. Indeed, 
lYusef-ZadeheFaLl (|2()11) find QPOs in Sgr A* VLA Hght 
curve at 43 GHz with the amplitude AF ^ 0.03 Jy, which is 
similar to our predicted AF ^ 0.05 Jy at 87 GHz. Yet, oscilla- 
tions are prominent on only 1 of 5 nights in VLA data, which 
may indicate intermittency of QPOs in the Galactic Center. 
Our modeling provides a potential reason for such intermit- 
tency: QPOs are practically absent when the magnetic field is 
weak due destruction by magnetic field reversals. 

The total flux QPO amplitude AF w 0.15 Jy is larger than 
the LP flux amplitude AF ?» 0.06 Jy. However, the obser- 
vational error of the total flux can also be larger SMA can 
measure the LP flux up to the leakage level of ~ 0.3%, while 
the total flu x is only measured to ~ 1% due to calibration un- 
certainties (IMarrone et al.ll2007l) . Then it is about equally dif- 
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t[M] 

Fig. 5. — Accretion rate dependence on time M(t) (top) and normalized 
spectrogram of the total flux light curve (bottom) for simulation intervals 
with Af = 600M. The normalized spectrogram is made up of ratios of 
log-smoothed to O.OSdex periodograms over log-smoothed to 3.0dex peri- 
odograms. The higher ratio and the lighter color indicate QPOs. 



ficult to detect total flux oscillations and LP flux oscillations. 
ALMA gives more hope than older observatories in detect- 
ing Sgr A* sub-mm QPOs. ALMA covers a wide range of 
sub-mm frequen cies 84 — 720 GHz and has a collecting area 
of -- 7 X 10^^ m^ (Brown etal. 2004) about 30 times the col- 
lecting area of SMA. ALMA observations of Sgr A* are ex- 
pected to have a flux error under 0.05 Jy and sampling every 
few minutes, which should be enough to reveal the predicted 
oscillations were they present on an observation night. 
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